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INTRODUCTION 


Long  range  target  detection  via  ionospheric  propagation  modes  — 
commonly  known  as  over-the-horizon  (OTH)  detection  —  has  long  been  considered 
as  a  means  of  extending  radar  coverage  to  large  distances.  Many  experimental 
studies  have  been  conducted  over  the  years  and  at  the  present  time  operational 
systems  are  under  development.  This  mode  of  detection  possesses  several  dis¬ 
tinct  advantages  including  an  ability  to  detect  targets  to  distances  in  excess 
of  3000  km.  However,  it  also  suffers  from  shortcomings  which  have  tended  to 
retard  its  development.  Among  these  are  the  fact  that  ionospheric  propagation 
modes  require  the  frequency  of  operation  to  be  less  than  approximately  three 
times  the  ionospheric  critical  frequency  —  -  10  MHz  —  leading  to  a  somewhat 
complex  propagation  analysis  for  precise  target  location.  The  problem  is 
especially  severe  at  high  latitudes  where  the  ionosphere  is  highly-variable  as 
a  result  of  auroral  particle  precipitation.  In  addition,  particle  precipita¬ 
tion  may  produce  enhanced  ionization  of  the  D-region  of  the  ionosphere  which 
can  lead  to  severe  absorption  of  the  energy  in  ionospheric  propagation  modes 
and  an  associated  reduction  in  the  radar  sensitivity.  Finally,  radars  utili¬ 
zing  ionospheric  propagation  paths  are  especially  sensitive  to  ionospheric 
clutter  from  magnetic-field-aligned  electron  density  irregularities.  Radar 
signals  backscattered  from  these  irregularities  can  be  very  intense  and  any 
effort  to  mitigate  against  their  influence  requires  knowledge  of  at  least  the 
Doppler  shift  and  spread  that  they  introduce  in  the  backscattered  signals. 

In  this  report  we  concentrate  on  the  third  of  these  problems  associ¬ 
ated  with  OTH  detection.  Our  knowledge  of  the  Doppler  spectral  characteris¬ 
tics  of  HF  signals  backscattered  from  ionospheric  irregularities  has,  in 
general,  been  quite  poor.  This  has  been  especially  true  for  backscatter  from 
F-region  irregularities  at  high  latitudes.  While  experimental  studies  at  high 


latitudes  have  been  conducted  the  equipment  has  often  been  designed  for  OTH 
detection  and  not  for  an  analysis  of  the  sp. _tral  characteristics  of 
ionospheric  irregularities.  Consequently,  when  Elkins  (1980)  produced  a  model 
for  HF  radar  auroral  clutter,  he  was  forced  to  assume  certain  spectral 
properties  for  the  signals  backscattered  from  high  latitude  irregularities. 
We  shall*  show  in  this  report  that  these  assumptions  were  of  limited 
validity.  In  particular,  we  shall  show  that  the  variability  in  the  magnitude 
and  direction  of  the  irregularity  drift  as  well  as  the  variability  in  spread 
of  the  Doppler  returns  make  it  difficult  to  discriminate  against  ionospheric 
clucter  by  Doppler  techniques,  alone. 

The  data  presented  in  this  report  have  been  obtained  with  The  Johns 
Hopkins  University/Applied  Physics  Laboratory  (APL)  HF  Radar  located  at  the 
Air  Force  Geophysics  Laboratory  High  Latitude  Ionospheric  Observatory  in  Goose 
Bay,  Labrador.  This  is  an  auroral  zone  site  and  signals  from  the  radar  are 
subjected  to  all  of  the  detrimental  ionospheric  effects  described  in  the  pre¬ 
vious  paragraphs.  The  radar  was  initially  funded  by  the  Atmospheric  Sciences 
Division  of  the  National  Science  Foundation  and  by  the  Atmospheric  and  Chemi¬ 
cal  Sciences  Directorate  of  Che  Air  Force  Office  of  Scientific  Research  for 
Che  study  of  small-scale  ionospheric  irregularity  structure.  Subsequently, 
additional  support,  including  a  second  receiving  array  for  the  antenna  system, 
have  been  provided  by  the  Electromagnetic  Propagation  Branch  of  Rome  Air 
Development  Center  and  by  the  Ionospheric  Effects  Division  of  the  Defense 
Nuclear  Agency.  At  the  present  time  the  radar  is  operating  at  a  modest  power 
level  Chat  is  approximately  30  dB  below  that  required  for  OTH  detection. 
Improvements  presently  being  made  to  the  radar  will  increase  the  sensitivity 
by  approximately  10  dB. 
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Since  the  APL  radar  was  developed  for  the  investigation  of  small- 
scale  ionospheric  structure,  it  has  several  capabilities  that  differentiate  it 
from  radar  systems  that  have  been  utilized  for  OTU  detection.  Most  notable  of 
these  is  an  unusual  multipulse  transmission  pattern  that  yields  unambiguous 
determinations  of  the  autocorrelation  functions  of  the  backscattered  signals 
as  a  function  of  range.  The  autocorrelation  functions  are  calculated  in  real 
time  by  an  on-line  microcomputer  and  stored  on  digital  tape.  Subsequent 
processing  at  APL  yields  Doppler  spectra  of  the  backscattered  signals.  The 
approach  as  it  is  implemented  on  the  radar  provides  for  a  Doppler  bandwidth  of 
330  Hz  about  the  transmitted  frequency.  At  an  operating  frequency  of  10  MHz, 
this  bandwidth  yields  unambiguous  line-of-sight  velocity  measurements  up  to 
nearly  2500  m/s. 

In  the  following  we  present,  first,  a  discussion  of  various 
procedures  that  one  might  use  for  spectral  and/or  autocorrelation  analysis  of 
backscattered  radar  signals.  This  discussion  is  derived  from  a  similar 
analysis  presented  in  an  earlier  RADC  report  (Greenwald,  1982).  Next,  we 
describe  the  method  by  which  the  autocorrelation  functions  obtained  by  the 
radar  are  processed  at  APL  in  order  to  remove  contributions  from  ground 
backscatter  and  unwanted  signal  sources.  The  autocorrelation  functions  are 
then  analyzed  to  determine  their  decorrelation  time  (inverse  of  the  spectral 
bandwidth)  and  their  period  (inverse  of  the  Doppler  frequency).  Finally,  the 
parameters  are  converted  to  Doppler  velocity  width  and  shift  through  inversion 
and  correction  for  frequency  of  operation.  These  values  have  been  used  in  a 
statistical  study  of  the  occurrence  and  Doppler  characteristics  of  backscatter 
from  high  latitude  ionospheric  irregularities.  The  results  from  this  analysis 
are  then  discussed  in  terms  of  their  implications  on  the  operation  of  OTH 
radars  in  high  latitude  environments. 
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2.0  TECHNIQUES  FOR  DOPPLER  ANALYSIS  OF  IONOSPHERIC  CLUTTER 

Small-scale  ionospheric  structure  falls  within  a  class  of  scatterer 
that  is  known  as  a  distributed  soft  target.  This  type  of  scatterer  often  has 

a  rather  small  volume  scattering  cross-section,  often  10“*®  -  10“*^  m^/m^; 

1 3 

however,  as  the  total  scattering  volume  may  be  quite  large,  for  example  10 

m^,  the  signals  backscattered  from  the  volume  may  be  appreciable.  The  volume 

scattering  cross-section  associated  with  soft  targets  is  determined  from 

Fourier  decomposition  of  the  medium.  Since  the  radar  signals  undergo  Bragg 

scattering  by  the  irregularities,  the  critical  Fourier  component  is  given  by 

k  *  k-  „  -k  _ „ . .  In  the  case  of  a  backscatter  radar,  k- _  and  k  „„„  are  equal 

inc  scat  me  scat  ^ 

and  oppositely  directed;  hence  the  Fourier  component  of  interest  has  a  wave¬ 
length  that  is  equal  to  half  the  radar  wavelength.  The  amplitude  of  this 
critical  Fourier  component  is  determined  by  the  physical  processes  occurring 
in  the  medium.  Examples  of  how  it  is  obtained  are  given  by  Farley  (1972)  and 
in  the  case  of  HF  backscatter  by  Walker  et  al.  (1986).  A  detailed  discussion 
of  irregularity  scattering  cross-sections  is  beyond  the  scope  of  this  report. 

Let  us  now  consider  the  phase  variations  of  a  signal  backscattered 
from  a  given  volume  of  ionospheric  irregularities.  If  one  compares  the  back- 
scattered  signal  with  a  reference  oscillator  operating  at  the  frequency  of  the 
transmitted  signal,  the  phase  difference  between  the  two  signals  will  vary 
with  time.  This  phase  variation  is  the  Doppler  effect  due  to  the  bulk  motion 
of  the  scattering  volume.  In  addition,  one  will  note  that  the  difference 
signal  randomizes  or  loses  phase  memory  after  a  relatively  short  time.  This 
time  is  the  decorrelation  time  of  the  medium  and  it  transforms  to  the  width  of 
the  Doppler  spectrum  (see  Section  3.2).  Precisely  stated,  the  decorrelation 
time  is  related  to  the  decay  of  the  two-point  time-dependent  cross-correlation 
function  of  the  process  producing  the  irregularity  structures.  Physically,  it 
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may  be  thought  of  as  due  to  the  finite  lifetimes  of  the  individual 
irregularity  structures  and  the  random  or  turbulent  velocities  that  exist  in 
the  scattering  volume. 

The  question  arises  as  to  how  one  might  analyze  the  signals  returned 
from  the  scattering  volume.  In  this  report  we  consider  four  different 
approaches  using  pulsed  transmissions  (These  as  well  as  frn-cw  transmissions 
have  been  discussed  in  an  earlier  report  by  Greenwald  (1982)).  The  approaches 
have  been  selected  so  as  to  demonstrate  that  increasingly  conplex  procedures 
will  enable  one  to  obtain  unambiguous  information  from  increasingly  complex 
situations.  The  four  approaches  are: 

1)  long  pulse, 

2)  repetitive  pulse, 

3)  variable  lag  double  pulse,  and 

4)  multiple  pulse. 

We  consider  each  of  these  pulse  modes  as  they  are  implemented  for  HF  radar 
systems.  We  shall  assume  that  ionospheric  backscatter  may  occur  simultan¬ 
eously  over  a  range  interval  of  1500  km  (10  ms  in  radar  group  delay),  that  it 
can  exhibit  Doppler  shifts  in  excess  of  1200  m/s  (100  Hz  at  12.5  MHz)  and  that 
it  has  a  typical  decorrelation  time  of  30  ms. 

2.1  Long  Pulse 

Doppler  analysis  of  the  signals  scattered  from  a  single  long  trans¬ 
mitted  pulse  is  performed  as  shown  in  Figure  la.  For  this  transmission  mode 
the  length  of  the  transmitted  pulse  should  approximate  the  decorrelation  time 
of  the  medium.  Hence  we  assume  it  to  have  a  width  of  30  ms.  The  signals 
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scattered  by  the  ionospheric  structure  are  assumed  to  be  referenced  to  the 

transmission  frequency  by  passing  them  through  a  phase  coherent  receiver.  The 

» 

quadrature  outputs  of  the  receiver  contain  the  desired  Doppler  information. 
In  order  to  obtain  the  desired  spectral  bandwidth  the  receiver  must  have  a 
bandwidth  of  at  least  100  Hz  and  its  outputs  must  be  sampled  at  twice  that 
rate.  The  sampled  data  forms  a  time  series  which  may  be  denoted  as 

C(j)  =■  A(j)  ♦  iB( j)  1  <  j  <  N  (1) 

where  A(j)  and  B(j)  represent  the  jth  sample  of  the  quadrature  outputs  of  the 
receiver.  This  time  series  may  be  Fourier  analyzed  or  it  may  be  analyzed  to 
yield  the  autocorrelation  function  of  the  backscattered  signal.  We  shall 
determine  the  latter  quantity  to  maintain  uniformity  between  the  various 
analysis  procedures.  It  is  given  by 

R(k)  *  <C(j)C*( j+k)>  1  <  k  <  N-l  (2) 

where  the  asterisk  represents  the  complex  conjugate.  The  Doppler  spectrum 
associated  with  R(k)  may  be  obtained  by  Fourier  transformation. 

The  disadvantage  of  the  long  pulse  approach  should  be  quite  evident. 
Due  to  the  length  of  the  transmission,  backscattered  signals  are  received 
simultaneously  from  the  entire  1500  km  extent  of  the  scattering  region.  Con¬ 
sequently,  there  is  no  way  to  resolve  any  range-dependent  Doppler  velocity 
structure  that  may  be  present. 
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Repetitive  Pulse 


In  order  to  retain  range-dependent  Doppler  information  a  somewhat 
more  sophisticated  pulse  transmission  technique  is  required.  For  example,  it 
can  be  noted  that  the  long  pulse  is  transmitted  for  the  entire  length  of  the 
sampling  sequence.  This  is  not  necessary  and,  in  fact,  one  may  replace  the 
long  transmission  by  a  sequence  of  short  transmitted  pulses  with  one  pulse 
preceding  each  sample  of  the  received  signal  (See  Figure  lb).  With  this  pulse 
scheme,  the  bandwidth  of  the  receiver  is  matched  to  the  length  of  the 
transmitted  pulse  and  the  spatial  resolution  of  the  measurement  is  determined 
by  the  latter  quantity.  If  the  transmitted  pulse  were  to  have  a  typical 
length  of  100  ps,  then  the  range  resolution  would  be  15  km. 

With  the  repetitive  pulse  technique  the  sampled  time  series  and  the 
resulting  autocorrelation  function  may  be  written  as 

C(t,j)  -  A(t,j)  +  iB ( t , j )  1  <  j  <  N  (3) 


and 


R(t,k)  -  <C(t,j )C*(t, j+k»  1  <  k  <  N-l  (4) 

respectively.  Here,  t  represents  the  time  delay  between  transmission  of  a 
pulse  and  sampling  of  the  backscattered  signal.  Since  the  transmitted  pulse 
length  is  short  in  comparison  to  the  time  -between  pulses  (5  ms  for  a  200  Hz 
sampling  frequency),  it  is  entirely  possible  to  analyze  the  backscattered 
signals  from  many  different  ranges,  simultaneously. 

Unfortunately,  the  repetitive  pulse  technique  has  a  significant 
shortcoming.  Initially,  it  had  been  assumed  that  the  scattering  region  had  a 
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range  extent  of  1500  km.  For  the  present  discussion  we  will  assume  it  to 
extend  from  1000  to  2500  km.  If  t  is  set  so  that  the  data  is  sampled  at  a 
range  of  1000  km  (t  ■  6.7  ms),  then  the  sample  will  also  be  of  signals  back- 
scattered  from  the  immediately  preceding  pulse  due  to  ionospheric  structure  at 
a  delay  of  11.7  ms  and  a  range  of  1750  km.  In  fact,  for  every  delay  t  that 
one  might  select  there  are  two  distinct  ranges  that  can  contribute  to  the 
backscattered  signal.  In  practice  the  problem  may  be  even  more  severe  than 
demonstrated  here  since  sampling  frequencies  in  excess  of  200  Hz  may  be 
required  and  scattering  regions  extended  for  more  than  1500  km  in  range  are 
entirely  possible. 

A  particularly  good  example  of  the  repetitive  pulse  technique  as  it 
is  implemented  at  VHF  frequencies  for  studies  of  the  E-region  radar  aurora  is 
given  by  Balsley  and  Ecklund  (1972). 


2.3  Variable-Lag  Double  Pulse 

In  order  to  avoid  the  problem  of  range  ambiguity,  while  maintaining 
the  requisite  spectral  bandwidth  and  resolution,  it  is  necessary  to  utilize 
some  non-repetitive  multiple  pulse  technique.  The  simplest  of  these  is  the 
variable  lag  double  pulse  method  (See  Figure  lc).  In  this  approach  N  pairs  of 
pulses  are  transmitted  with  a  spacing  between  pulses  that  is  some  multiple  k 
of  an  elemental  spacing  T.  The  time  between  pulse  pairs  must  be  sufficiently 
long  that  no  signals  are  received  from  the  previous  pulse  pair  while  the  data 
from  the  current  pulse  pair  is  being  sampled.  By  referencing  the  time  delay  t 
to  the  time  of  the  first  pulse  of  a  pulse  pair,  we  can  write  the  sampled 
returns  from  the  transmissions  as 


C(t,j)  -  A( t , j )  +  iB(t,J) 
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and 


C(t+kT,j)  =  A(t+kT,j)  +  iB(t+kT,j) 


(6) 


If  we  further  allow  the  index  k  to  be  stepped  through  a  range  of  values  from  0 
to  M,  we  can  write  the  autocorrelation  function  derived  from  this  analysis 
procedure  as 


R(t ,k)  *  <C(t, j)C*(t+kT,j)>  0  <  k  <  M  (7) 


Let  us  consider  the  application  of  the  double  pulse  approach  to 
backscatter  from  extended  regions  of  ionospheric  structure.  If  we  denote  the 
subscript  1  as  indicating  a  response  due  to  the  first  transmitted  pulse  and 
the  subscript  2  as  due  to  the  second  transmitted  pulse,  7)  can  be  rewritten  as 

R(t,k)  *  <(CX  (t,j)+C2  (t,j))(C1*  (t+kT, j)+C2*  (t+kT , j ) )>  (8) 

Now,  by  referencing  the  response  of  each  transmitted  pulse  to  the  time  of  that 
pulse,  explicitly  noting  the  time  delay  between  the  two  samples,  and  expanding 
the  expectation  value,  (8)  may  be  rewritten  as 

R(t,k)  *  <CX  (t.j.OC^  (t+kT, j,k)»<C1  (t,j,0)C2*  (t,j,k)> 

+  <C2  (t-kT, j,0>Cl*  (t+kT, j,k)>+<C2  (t-kT, j,0)C2*  (t,j,k)> 

0  S  k  <  M  (9) 
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Examination  of  (9)  shows  that  only  the  second  term  on  the  right-hand-side 
involves  a  comparison  of  data  from  the  same  group  delay.  All  other  terms 
average  to  zero  since  the  returns  from  two  distinct  ranges  are  uncorrelated. 
They  contribute  to  the  autocorrelation  function  only  by  raising  the  overall 
noise  background  of  the  measurement. 

As  in  the  case  of  the  repetitive  pulse,  the  variable  lag  double 
pulse  method  may  be  used  for  simultaneous  examination  of  the  Doppler  charac¬ 
teristics  of  signals  backscattered  from  many  different  ranges.  Moveover,  the 
elemental  lag  T  and  the  number  of  lags  M  may  be  adjusted  so  that  one  can 
always  perform  a  complete  and  correct  Doppler  analysis  on  the  backscattered 
signals.  Unfortunately,  the  method  is  very  inefficient  in  terms  of  processing 
time.  For  the  longer  lags  of  the  autocorrelation  function  one  must  wait  for 
many  tens  of  milliseconds  to  obtain  the  appropriate  pulse  separations.  Since 
the  lags  are  performed  sequentially  and  N  must  be  sufficiently  large  to  obtain 
good  statistics,  it  is  highly  questionable  whether  the  high  latitude  iono¬ 
sphere  is  sufficiently  time  stationary  for  this  method  to  be  used  in  the  HF 
frequency  regime.  An  example  of  its  application  to  radar  auroral  backscatter 
at  upper  VHF  frequencies  is  given  by  Nielsen  et  ai.(1984). 

2.4  Muitipulse 

It  would  be  highly  desirable  to  maintain  the  positive  features  of 
the  double  pulse  technique  while  improving  its  efficiency.  One  might  ask 
whether  it  is  possible  to  transmit  multiple  pairs  of  double  pulses  with 
different  lags  in  such  a  manner  that  the  pulse  pairs  do  not  interfere  with  one 
another.  This  question  has  been  answered  in  the  affirmative  and  the  solutions 
form  the  basis  of  muitipulse  transmission  techniques  (See  Farley,  1972).  With 
muitipulse  techniques  one  seeks  to  transmit  a  pattern  of  pulses  such  that  the 
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pulse  separations  contain  most  or  all  of  the  lags  necessary  for  determining 
the  autocorrelation  function  of  the  backscattered  signal  and  such  that  there 
is  little  or  no  ambiguity  in  the  ranges  that  contribute  to  any  given  lag.  In 
practice,  a  four  pulse  transmission  pattern  is  the  longest  that  can  be  used  to 
obtain  a  full  seven  lag  (counting  the  zero  lag)  autocorrelation  function  with 
no  range  ambiguity.  Patterns  with  a  greater  number  of  pulses  either  have  gaps 
in  their  associated  autocorrelation  patterns  or  they  have  range  ambiguity  for 
certain  lags. 

Farley  (1972)  has  suggested  several  six-pulse  patterns  that  might  be 
used  to  determine  fifteen  lags  of  a  17-lag  (counting  the  zero  lag)  autocorre¬ 
lation  function.  For  the  Goose  Bay  radar,  we  have  adopted  a  7-pulse  pattern, 
shown  in  Figure  Id,  that  is  used  to  determine  a  full  17-lag  autocorrelation 
function.  This  pattern  exhibits  range  ambiguity  for  k  =  1,2, 3, and  13; 
however,  the  large  time  separation  between  the  two  repetitions  of  k  *  1,2, and 
3,  makes  it  extremely  unlikely  for  two  distinct  range  intervals  to  yield 
simultaneous  backscattered  signals.  The  remaining  ambiguous  lag  is  a 
problem.  Either  it  can  be  dropped  from  the  analysis  or  it  can  be  ignored  due 
to  the  small  residual  correlation  that  exists  in  the  cross  product  analysis  at 
these  large  lags. 

In  addition  to  retaining  the  benefits  of  the  variable  lag  double 
pulse  analysis,  the  multipulse  technique  is  quite  efficient  from  an  analysis 
point  of  view.  With  it,  the  data  for  a  full  17-lag  autocorrelation  function 
may  be  acquired  in  16T.  In  contrast,  at,  least  300T  would  be  required  to 
acquire  the  equivilent  data  using  the  variable  lag  double  pulse  technique. 
Thus,  there  is  a  speed  enhancement  of  more  than  a  factor  of  18  in  the  data 
acquisition. 


11  - 


Some  minor  shortcomings  of  the  multipulse  approach  are  an  increased 
noise  background  due  to  uncorrelated  scatter  from  unwanted  ranges  and  the 
occurrence  of  transmitter  pulses  at  certain  delays  which  affect  the  autocor¬ 
relation  analysis.  Fortunately,  neither  of  these  problems  is  particularly 
serious. 


3.0  PROCESSING  OF  AUTOCORRELATION  FUNCTIONS  FROM  THE  GOOSE  BAY  RADAR 

In  the  case  of  the  APL  Goose  Bay  radar,  the  analysis  of  the  multi¬ 
pulse  autocorrelation  functions  is  carried  out  by  a  FORTRAN  program,  FITACF. 
This  is  a  rather  complex  program,  but  its  primary  function  is  to  fit  a 
theoretical  autocorrelation  function  to  the  observed  data  in  order  to  deter¬ 
mine  the  true  backscattered  power,  the  Doppler  velocity,  and  the  spectral 
width.  Each  5  second  integration  is  separately  processed  in  two  basic  steps. 
Noise  Reduction  and  Parameter  fitting  of  each  individual  ACF. 

3.1  Noise  Reduction 

There  are  several  sources  of  noise  which  complicate  the  analysis  of 
the  radar  data.  The  cosmic  HF  background  noise  is  relatively  constant  in  time 
but  has  a  frequency  dependence  of  the  form  P  -  f^  5/2) ^  The  muicipuise  tech¬ 
nique,  as  noted  above,  causes  strong  scatter  at  some  ranges  to  contribute  to 
the  noise  level  at  other  ranges.  In  addition,  there  are  local  sources  of  HF 
noise  produced  by  nearby  equipment,  as  well  as  the  inherent  noise  of  the 
receiver  and  digitizers.  Finally,  in  the  HF  frequency  band,  there  are  nearly 
always  remote  radio  transmitters  that  contribute  to  the  noise  background  at 
varying  intensity  levels.  Some  of  these  transmitters  are  nearly  CW  sources 
while  others  are  highly  modulated  or  pulsed. 
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The  first  stage  of  the  noise  reduction  process  is  the  determination 


of  the  basic  noise  level  and  the  "noise  ACF."  An  initial  noise  level  is 
determined  from  the  average  backscattered  lag-0  power  from  the  10  weakest 
ranges.  An  average  noise  ACF  is  then  formed  from  all  the  autocorrelation 
functions  which  have  lag-0  power  less  than  1  dB  above  the  initial  estimate  of 
the  noise  level.  In  a  typical  case  where  there  are  no  external  transmitters 
producing  noise,  the  noise  ACF  will  have  a  non-zero  power  at  lag-0  and  be 
nearly  zero  for  all  other  lags.  If  a  CW  transmitter  is  present,  however,  it 
will  be  present  at  all  range  gates.  In  this  case,  the  noise  ACF  will  show  a 
clear  non-zero  frequency.  In  either  case,  the  noise  ACF  is  then  subtracted 
from  the  raw  ACFs.  This  subtraction  removes  the  excess  lag-0  power  due  to 
pure  random  noise  and  may  also  remove  or  substantially  reduce  the  effect  of 
other  coherent  noise  sources  such  as  CW  transmitters. 

To  see  the  effect  of  removing  the  noise  ACF  more  clearly,  let  us 
consider  what  the  combination  of  two  separate  sources  produces  for  an  auto¬ 
correlation  function.  For  simplicity  we  shall  assume  each  source  is  perfectly 
correlated  with  itself,  but  is  not  correlated  with  the  other  source.  This 
will  be  the  case  when  one  source  is  an  external  transmitter  and  the  other  is  a 
backscattered  ionospheric  signal  with  a  very  narrow  Doppler  spectrum.  In  such 
a  case  the  autocorrelation  function  produced  by  a  single  transmission  of  the 
multipulse  sequence  will  be 
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where  $  is  an  arbitrary  phase.  When  many  pulse  sequences  are  added  up,  the 
cross-term  will  tend  to  die  out  since  the  phase  is  different  each  time.  The 
expected  result  is 
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A2  e  1  ♦  B2  e  2  ♦  J£-  ,  1  2  e1*  (11) 

where  N  is  Che  number  of  pulse  sequences  produced  during  Che  5  second  integra¬ 
tion.  The  value  of  H  is  typically  around  60  and  hence  the  cross  term  is 
suppressed  by  a  factor  of  7.7.  Let  us  assume  that  the  kl  term  is  a  noise 
transmitter  and  the  8*  term  is  the  ionospheric  signal.  If  A  >  B  then  on 
subtracting  the  A  term  (the  noise  ACF)  from  the  total  ACF  we  are  left  with 


82  e  2  ♦  -£L  e  1  2  e1* 


If  we  are  to  be  able  to  extract  the  true  ionospheric  ACF  from  the  remaining 
noise  we  must  have 


In  terms  of  power  this  may  be  rewritten  as 

<  N  (14) 

B 

For  our  case,  where  N  *  60  this  implies  that  the  ionospheric  signal  can  be 
detected  and  extracted  from  the  noise  even  when  it  is  17  dB  below  the  noise. 
Of  course  this  is  an  ideal  case  and  still  requires  an  ionospheric  signal 
larger  than  the  true  random  noise  present  which  in  actual  practice  limits  the 
usefulness  of  this  technique  to  cases  where  the  noise  source  and  the  iono¬ 
spheric  signal  are  comparable. 
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3.2  Parameter  Fitting 

Once  the  noise  ACF  has  been  removed  from  the  data,  it  is  possible  to 
fit  an  assumed  functional  form  to  each  of  the  ACFs.  We  have  looked  at  two 
possible  functional  forms,  an  exponential  decorrelation  with  time  and  a 
Gaussian  decorrelation.  If  we  assume  an  exponential  decorrelation  we  will 


_  _  luit  -Xt 

V‘>  *  cx  *  * 


where  is  the  power,  u  is  the  Doppler  frequency  and  X  is  the  decorrelation 
parameter.  The  Fourier  transform  of  gives  a  Doppler  spectrum  of  the  form 

V“'>'  72—7-77  (16) 

X  +  (a»-w  ) 

which  peaks  at  u"  a  u  and  has  a  full  width  at  half-maximum  of  2X  radians/s. 
Alternatively,  if  we  assume  a  Gaussian  decorrelation  we  will  have 


R  (t)  =  C  eiwt  e-*7^2 
a  o 


and  the  Doppler  spectrum  will  be 


-(u-u' 
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S  (u  )  =  —  C  e 
0  00 


which  again  peaks  at  m'  -  to  and  has  a  width  of  4  o  /  In  (2)  5  2.76  c.  The 


parameters  C^,  X,  u  or  C^,  a,  u  are  fitted  using  a  power-weighted  least- 
squares  fit.  Weighting  by  the  power  emphasizes  the  importance  of  the  earlier 
lags  where  the  power  is  greater  and  the  errors  and  influence  of  noise  are 


less.  The  Doppler  frequency  is  found  by  doing  a  least-squares  fit  to  the 
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observed  phase  at  each  lag  (Hanuise  et  al.,  1985),  while  Che  decorrelation 
parameter  and  power  (C^  and  X  or  C^and  a)  are  found  by  doing  a  fit  to  the 
logarithm  of  the  magnitude  of  the  observed  ACF. 

There  are  several  complications  in  doing  the  fit  that  must  be  con¬ 
sidered.  The  primary  difficulty  is  the  problem  of  "bad  lags".  Since  some 
samples  are  taken  at  the  same  time  that  the  transmitter  is  on  and  the  receiver 
is  off,  a  range  which  uses  one  of  these  samples  in  computing  its  ACF  will  have 
one  or  more  lags  which  are  bad  and  should  not  be  used  in  the  fitting  process. 
In  addition  to  lags  which  are  bad  because  of  transmitter  pulses,  lag-13  has 
multiple  redundancy  in  the  multipulse  sequence  and  can  lead  to  range-aliased 
results.  Finally,  pulsed  signals  from  remote  signal  sources  can  cause  sudden 
changes  in  the  received  signal  which  affect  only  a  few  samples.  The  noise 
from  these  pulsed  sources  cannot  be  removed  by  the  method  described  in  the 
previous  section,  and  each  lag  which  uses  one  of  the  affected  samples  will  be 
a  bad  lag. 

The  first  step  in  performing  the  fit  to  the  observed  data  is  there¬ 
fore  to  detect  bad  lags  and  eliminate  them  from  the  data  used  for  the  fit. 
The  samples  which  are  taken  during  transmitter  pulses  are  readily  determined 
and  removed  from  the  autocorrelation  function  analysis.  In  addition,  because 
of  the  range  aliasing  problem,  lag-13  for  all  ranges  is  considered  a  bad 
lag.  The  remaining  lags  are  scanned  and  any  lag  which  shows  a  large  increase 
in  magnitude  over  the  preceding  lags  is  considered  a  bad  lag. 

Once  the  bad  lags  have  been  removed  the  remaining  lags  are  used  to 
determine  the  parameter  fit.  Although  the  determination  of  C  and  X  (or  C 
and  o)  is  straightforward,  there  is  one  additional  complication  in  determining 
the  Doppler  frequency.  To  do  a  straight  line  fit  to  the  phase  we  assume  that 
the  phase  is  either  monotonically  increasing  or  decreasing  (or  approaching 
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receding  Doppler  velocities).  The  actual  measured  phase ,  however,  is  confined 
to  the  range  ±  it.  Multiples  of  2it  must  therefore  be  added  to  the  observed 
phase  at  the  higher  lags.  The  basic  approach  used  is  described  by  Hanuise,  et 
al.  (1985),  but  it  has  been  modified  to  improve  the  determination  of  the 
2ir  folding  factor  in  cases  where  several  lags  are  bad  and  cannot  be  used.  In 
addition,  the  error  in  determining  id  is  estimated  and  used  to.  determine  the 
velocity  error. 

The  question  of  which  of  the  two  functional  forms  for  the  ACF  (Equa¬ 
tions  15  and  17)  is  the  more  appropriate  one  to  use  is  difficult  to  determine. 
If  the  decorrelation  is  due  to  the  fact  that  we  are  observing  multiple  scat- 
terers  in  the  presence  of  velocity  turbulence,  then  we  would  expect  the 
Doppler  spectrum  to  be  Gaussian  and  Equations  (17)  and  (18)  to  be  the  more 
appropriate  ones.  If,  on  the  other  hand,  the  decorrelation  is  due  to  the 
growth  and  decay  of  the  plasma  irregularity  structures  with  little  or  no 
velocity  turbulence,  then  the  appropriate  equations  are  (15)  and  (16). 

To  investigate  this  further,  we  have  tried  both  fits  on  several 
different  periods  of  data.  An  example  is  shown  in  Figure  2.  The  data  were 
taken  on  January  6,1986  at  20:01:40  UT  (Figure  2a)  and  20:01:45  UT  (Figure  2b) 
for  the  same  range  (1080  km).  Table  1  gives  the  values  of  the  parameter  fit 
for  the  two  methods  and  the  power  weighted  error  for  each. 
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TABLE  1 


log(Cx) 

X 

CX 

width  (  X) 


20:01:40 

11.71 

48.89 

0.042 

97.8  rad/s 


20:01:45 

11.51 

43.11 

0.022 

86.2  r ad/s 


log(Co) 

o 

c 

0 

width  (a) 


11.48 

38.56 

0.007 

106.81  rad/s 


11.25 
31.51 
0.039 
87.3  rad/s 


For  the  first  case  (Figure  2a),  the  Caussian  fit  is  clearly  superior,  but  the 
physical  results,  power  and  width,  are  not  markedly  different.  In  the  second 
case  (Figure  2b),  the  exponential  fit  has  the  smaller  error,  but  again,  the 
physical  results  are  very  similar.  Since  these  two  examples  are  from  the  same 
range  and  are  separated  by  only  five  seconds,  it  is  unlikely  that  the  physical 
process  responsible  for  the  decorrelation  could  have  changed  greatly.  We 
conclude  that  from  the  data  currently  available,  it  is  nearly  impossible  to 
determine  which  decorrelation  model  is  the  "correct"  one  and  the  power  and 
spectral  width  determined  by  either  method  are  very  similar.  For  this  statis¬ 
tical  survey,  we  have  used  the  exponential  fit  represented  by  Equation  (15). 


3.3  Separation  of  Ground  Scatter  and  Ionospheric  Scatter 

Once  the  parameter  fit  has  been  performed  on  an  ACF,  it  is  possible 
to  separate  ground  scatter  from  ionospheric  scatter.  Ground  scatter  is  char¬ 
acterized  by  a  very  low  Doppler  frequency  and  a  very  small  decorrelation 
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parameter  (i.e.  large  decorrelation  time).  After  analyzing  several  hundred 
ACFs  from  many  different  periods  we  have  determined  a  reasonable  set  of  limits 
which  define  the  spectra  of  ground  backscattered  signals: 

1)  Ground  backscattered  signals  typically  have  |u)j  <  18.85 

radians /s  (3  Hz). 

2)  Ground  backscattered  signals  have  small  |x  | .  The  actual  limit 
depends  on  the  power  (C).  As  can  be  seen  in  Figure  3,  weak 
ground  scatter  signals  tend  to  have  larger  values  of  X  than 
strong  signals.  The  limits  we  have  used  here  are 

5  s”1  for  SNR  >  6  dB 
(12-SNR)  s"1  for  0  <  SNR  <  6  dB 
12  s-1  for  SNR  <  0  dB 

With  these  choices  for  the  limits  defining  ground  scatter  we  find  that  we 
misidentify  approximately  5-102  of  weak  ground  scatter  signals  as  being  true 
ionospheric  signals.  Also,  a  small  percentage  of  weak  ionospheric  signals, 
which  have  very  low  velocities  are  misidentif ied  as  ground  scatter,  but  it  is 
difficult  to  determine  exactly  what  this  percentage  is.  Strong  signals  (>  10 
dB)  are  almost  never  misidentif ied. 

If  the  ACF  has  been  identified  as  due  to  ground  scatter,  the  ground 
scatter  function  defined  by  C,  u  and  X  is  subtracted  from  the  data  and  the 
residual  is  re-analyzed.  If  a  reasonable  fit  is  found  for  the  parameters  of 
the  residual  (i.  e.  |u|  >  18.85  and  relative  error  <  1)  the  ACF  is  identified 
as  a  mixture  of  ground  scatter  plus  ionospheric  scatter. 
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4.0 


STATISTICAL  METHOD 


4.1  Selection  of  the  Data 

In  order  to  insure  an  unbiased  sampling  of  the  radar  data,  all  the 
data  available  from  July  1,  1984  to  August  30,  1985  were  divided  into  1-hour 
sections.  A  random  sample  without  replacement  was  then  done  to  select  177  1- 
hour  periods  of  data.  From  July  1984  to  March  1985  the  radar  was  being  oper¬ 
ated  in  a  mode  which  allowed  us  to  collect  only  4  hours  of  data  per  day. 
During  some  of  this  7  month  period  the  4  hours  were  spread  evenly  over  the 
entire  day,  taking  only  10  minutes  of  data  each  hour.  At  other  times,  the 
data  were  taken  continuously  for  4  hours  and  the  radar  was  not  operated  for 
the  remainder  of  the  day.  Thus,  the  available  data  from  the  first  seven 
months  of  this  study  were  not  uniformly  distributed  and  this  must  be  taken 
into  account  when  we  analyze  the  statistics.  Starting  in  March  1985,  new 
software  was  installed  which  allowed  the  radar  to  operate  24  hours  a  day,  and 
the  data  from  the  final  6  months  of  this  study  are  nearly  uniformly 
distributed. 

Once  the  data  had  been  selected,  they  were  divided  into  two  cate¬ 
gories,  scatter  and  no-scatter.  A  1-hour  period  was  assigned  to  the  scatter 
category  provided  that  there  was  at  least  1  10-minute  period  within  the  data 
for  which  backscactered  power  was  clearly  present  with  a  signal-to-noise  ratio 
of  3  dB  or  more.  If  no  such  10-minute  period  could  be  found,  the  1-hour 
period  was  assigned  to  the  no-scatter  category.  The  10-minute  limit  was 
required  for  two  reasons;  firstly,  to  limit  the  total  amount  of  data  we  would 
have  to  process  and,  secondly,  because  there  were  many  1-hour  periods  when 
only  10  minutes  of  data  were  actually  taken.  It  should  be  stressed,  however, 
that  our  definition  of  the  probability  of  observing  scatter  within  a  given 
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hour  means  Che  probability  of  observing  at  least  10  minutes  of  backscattered 
power  at  or  above  3  dB  during  that  hour. 

Having  determined  which  of  the  1-hour  periods  show  good  backscat¬ 
tered  power,  we  then  selected  10  minutes  of  data  from  each  period  for  more 
detailed  analysis.  In  order  to  explain  the  techniques  used  in  the  remainder 
of  the  data  analysis,  it  is  first  necessary  to  describe  how  the  radar 
operates.  A  full  description  of  the  Goose  Bay  HF-Radar  can  be  found  in 
Greenwald,  et  al.  (1984).  The  radar  is  a  phased-array  system  which  is  elec¬ 
tronically  steerable  in  16  beam  directions.  Data  are  taken  for  a  single  beam 
for  an  integration  period  of  5  seconds.  A  full  scan  of  16  beams  thus  takes  80 
seconds.  There  are  50  non-overlapping  range-gates,  each  gate  being  30  km 
wide.  For  each  5  second  integration  the  autocorrelation  functions  are  ana¬ 
lyzed  by  FITACF .  The  combined  results  of  FITACF  for  all  the  data  may  then  be 
used  for  further  statistical  study. 

4 . 2  Event  Analysis 

The  data  from  FITACF  are  first  analyzed  in  terms  of  the  "events". 
By  "event"  we  mean  all  the  data  from  one  10-minute  selection  of  data.  The 
events  are  put  into  a  4-dimensional  contingency  table, which  is  defined  in 
Table  2.  Tables  of  reduced  dimension  are  produced  by  summing  over  one  or  more 
of  the  dimensions.  A  description  of  the  conputer  program  used  to  generate 
these  tables  is  given  In  the  appendix.  An  event  is  considered  to  be 
ionospheric  if  252  or  more  of  the  autocorrelation  functions  which  had  power 
greater  than  the  noise  level  (hereafter  simply  referred  to  as  ACFs)  fulfilled 
the  criteria  in  FITACF  to  be  considered  ionospheric  scatter.  Similarly,  an 
event  is  considered  to  be  ground  scatter  if  252  or  more  of  the  ACFs  fulfilled 
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Che  criteria  Co  be  considered  ground  scatter.  If  both  conditions  are  met  the 
event  is  categorized  as  "both".  A  similar  definition  is  used  in  determining 
Che  scaCter  region.  If  25Z  or  more  of  the  ACFs  come  from  ranges  between  300 
and  600  km  the  event  is  considered  to  be  E-region  scatter.  If  25Z  or  more  of 
the  ACFs  come  from  ranges  900  to  2100  km  the  event  is  considered  to  be  F- 

region,  if  both  conditions  are  met  the  event  is  described  as  "both"  and  if 

neicher  condition  is  met  (i.e.  the  ACFs  were  primarily  from  the  ranges  between 
600  and  900  km)  the  region  is  "undetermined".  Although  such  definitions  are 

meaningless  for  the  ground  scatter  events  the  use  of  the  "region"  category 

does  not  cause  a  problem  since  the  events  which  are  strictly  ground  scatter 
events  can  be  ignored  when  looking  at  lower  dimensional  contingency  tables. 
The  Kp  bins  are  in  integral  values  of  Kp  and  thus,  for  example,  bin  3  included 
Kp  of  3-,3,3+. 

TABLE  2 

Definition  of  the  "Event"  Table 


Quantity 

#  of  bins 

bin  ranges 

Time  (UT) 

4 

0-6,6-12,12-18,18-24 

Scatter  Type 

3 

ionospheric,  ground,  both 

Region 

4 

E,  F,  both,  undetermined 

Kp 

10 

0,1, 2, 3, 4, 5, 6, 7, 8, unknown 

After  examining  the  data  in  tepras  of  "events"  we  turn  to  the 
individual  ACFs.  First,  the  ACFs  are  divided  into  three  groups,  those  that 
have  ionospheric  backscattered  power  above  the  noise  level,  those  that  have 
ground  scatter  power  above  the  noise  Level,  and  those  which  show  no  power 
above  the  noise  level.  We  found  a  total  of  42,852  ACFs  which  were  ionospheric 
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As  described  in  Sections  4.1  and  4.2,  the  first  step  of  our  statist¬ 
ical  analysis  has  been  to  determine  the  probability  of  observing  appreciable 
backscattered  power.  The  selected  1-hour  periods  were  divided  into  scatter  or 
no-scatter  categories  and  were  binned  by  time  into  four  6-hour  intervals,  0  - 
6  UT,  6-12  UT,  12  -  18  UT,  and  18  -  24  UT.  The  magnetic  local  time  at  Goose 
Bay  is  approximately  3  hours  earlier  than  UT  and  these  four  time  bins  are 
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therefore  approximately  centered  around  magnetic  midnight,  dawn,  noon  and 
dusk.  The  results  are  summarized  in  Table  4  and  Figure  4.  The  rows  in  Table  4 
labeled  "obs."  give  the  actual  number  of  cases  observed  for  each  time  period, 
the  rows  labeled  "expect"  give  the  expectation  values  based  on  the  marginal 
probabilities.  The  sum  of  (expected-observed)*  /expected  is  the  chi-square 
statistic.  The  value  of  x  /4  (4  is  the  number  of  degrees  of  freedom) 
strongly  indicates  that  there  is  a  statistically  significant  time  dependence. 
Overall,  the  probability  of  observing  some  kind  of  scatter  (ionospheric  or 
ground)  is  approximately  582  with  the  highest  probability  (842)  occurring  in 
the  local  midnight  sector  and  the  smallest  (292)  in  the  local  noon  sector.  We 
next  turn  to  the  "event"  statistics.  In  order  to  improve  the  statistics  for 
the  event  table  and  ACF  table  we  randomly  selected  an  additional  23  periods  of 
data  that  had  some  sort  of  backscattered  power,  thus  bringing  the  total  number 
of  events  to  125.  The  first  table  to  consider  is  the  2-d  table  of  scatter 
type  vs  time.  Table  5  and  Figure  5  show  the  distribution  and  chi-square  test 
for  this  Cable.  Although  there  was  a  strong  time  dependence  on  the  overall 
probability  of  seeing  some  kind  of  scatter  (Fig.  4),  there  is  little 
dependence  of  the  type  of  scatter  on  time.  The  chi-square  test  shows  no 
significant  deviation  of  the  number  of  observed  events  from  the  expectation 
values. 

If  we  look  at  the  type  of  scatter  as  a  function  of  Kp,  however, 
quite  a  different  picture  emerges.  Table  6  and  Figure  6  show  the  distribution 
of  events  by  scatter  type  and  Kp.  Clearly, , ionospheric  scatter  as  opposed  to 
ground  scatter  increases  as  Kp  increases.  The  chi-square  test  confirms  the 
statistical  significance  of  this  result  with  over  992  confidence. 

The  next  question  is  how  the  ionospheric  scattering  region  is 
affected  by  time  and  Kp.  Table  7  and  Figure  7  show  the  distribution  of  the 


-  24  - 


events  by  region  and  time.  Note  that  the  events  where  the  region  could  not  be 
determined  have  been  ignored.  We  find  that  ionospheric  scatter  from  the  F- 
region  is  much  more  likely  in  the  time  range  from  12  to  24  UT  (i.e.  the  local 
noon  and  dusk  sectors).  Conversely,  E-region  scatter  is  significantly  more 
likely  during  the  dawn  (6-12  UT)  sector.  Table  8  and  Figure  8  show  the 
distribution  of  the  events  by  region  and  Kp.  Although  Figure  8  seems  to 
indicate  a  trend  of  increasing  probability  of  E-region  scatter  with  increasing 
Kp,  this  cannot  be  confirmed  by  the  chi-square  test.  However,  the  chi-square 
test  is  not  particularly  sensitive  to  trends  when  there  are  many  degrees  of 
freedom.  The  probability  of  short  wavelength  E-region  irregularities  being 
formed  is  known  to  be  related  to  the  strength  of  the  auroral  electrojets  which 
in  turn  are  related  (albeit  weakly)  to  Kp.  Thus,  although  the  chi-square  test 
cannot  confirm  the  trend  shown  in  Figure  7,  the  increase  in  the  percent  of  E- 
region  scatter  with  increasing  Kp  is  not  unexpected. 

4.4  Statistical  Distribution  of  Spectral  Observations 

We  turn  now  to  the  results  obtained  from  the  individual  spectra. 
First,  as  a  check  on  the  validity  of  Che  data  analysis  procedure  of  FITACF,  we 
can  compare  the  observed  lag-0  backscattered  power  with  the  value  obtained 
from  the  fit  to  the  full  autocorrelation  function.  When  the  backscattered 
power  is  due  to  ionospheric  scatter  these  two  quantities  should  be  highly 
correlated.  At  low  powers,  however,  we  would  expect  the  lag-0  power  to  be 
greater  Chan  the  fitted  power  since  the  lag-0  power  includes  some  residual 
noise  as  well  as  the  ionospheric  signal.  As  Figure  9  shows,  the  results  are 
as  expected.  The  points  used  are  all  the  ACFs  above  0  dB  which  were 
identified  as  ionospheric  scatter,  a  total  of  42,852  points.  The  analysis 
procedure  of  FITACF  includes  an  error  estimate  in  the  velocity  obtained  from 
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Che  lease-squares  fit  Co  Che  observed  phase  angle.  Figure  10  shows  Che 
variacion  of  Che  velocicy  error  as  a  funecion  of  Che  lag-0  power  (10a)  and  as 
a  funecion  of  Che  ficced  power  (10b).  We  would  expect  larger  velocicy  errors 
Co  occur  aC  lower  powers  and  this  is  indeed  whac  we  find  when  ploeeing  Che 
error  vs  Che  ficced  power.  However,  Che  scaCCer  ploc  obcained  by  ploCCing  Che 
error  vs  Che  lag-0  power  shows  a  clear  deviation  from  Che  expecced 
disCribucion.  Presumably,  chis  is  due  Co  noise  which  increases  boch  che  lag-0 
power  and  Che  velocicy  error.  The  consiscency  of  Che  resulcs  using  Che  ficced 
power  leads  us  co  conclude  chac  Che  ficced  power  is  Che  proper  paramecer  Co 
use  in  Che  remainder  of  our  scaciscical  analysis.  The  mean  value  of  che 
Doppler  velocity  error  for  che  daCa  shown  in  Figure  10b  is  Av  =  104  m/s,  while 
che  standard  deviation  is  a ^  -  141  m/s.  This  provides  a  reasonable  criterion 
for  rejecting  daCa  based  on  the  velocity  error.  A  velocity  error  greater  chan 
250  m/s  is  slightly  more  than  l-o  greater  than  the  mean  error  and  on  a  number 
of  our  scatter  plots  as  well  as  che  6-dimensional  contingency  table,  we  have 
restricted  the  data  Co  those  points  with  velocity  error  less  than  250  m/s.  It 
would  be  improper  Co  use  a  more  restrictive  limit  on  the  velocity  error 
because  che  error  is  correlated  to  the  Doppler  spectral  width.  Figure  11 
shows  a  scatter  plot  of  velocicy  error  vs  spectral  width.  It  is  clear  chat 
wide  spectra  tend  Co  give  rise  co  larger  errors  in  estimating  the  mean  Doppler 
velocity,  but  that  is  to  be  expected  and  the  wide  spectra  are  perfectly  valid 
data.  Thus,  using  too  stringent  a  velocity  error  criterion  for  accepting  data 
would  result  in  che  rejection  of  an  excessive  number  of  valid  data  having  wide 
spectral  characteristics. 

In  Figure  12  we  plot  the  spectral  width  as  a  function  of  fitted 
power.  The  relation  between  the  power  and  the  velocity  error  (Figure  10)  and 
the  relation  between  spectral  width  and  velocity  error  (Figure  11)  necessarily 
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implies  there  should  be  a  tendency  to  observe  larger  spectral  widths  at  lower 
powers.  Essentially,  the  increase  in  spectral  width  with  the  decrease  in 
backscattered  power  is  due  to  the  increasing  effect  of  noise  in  the  analysis 
of  the  autocorrelation  functions  at  lower  powers.  The  expected  inverse 
correlation  between  spectral  width  and  power  is  clear  in  Figure  12,  but  the 
distribution  shows  the  existence  of  larger  spectral  widths  at  all  powers  than 
the  velocity  error  distribution  in  Figure  10b.  The  partial  cut-off  in  the 
density  of  points  in  Figure  11  at  low  power  and  low  width  is  due  to  the 
classification  of  many  of  these  low  power  ACFs  as  coining  from  ground  scatter. 

We  next  turn  to  the  question  of  the  relation  between  ionospheric 
backscattered  power  and  range.  Walker,  et  al.  (1986)  have  shown  that  the 
intensity  of  the  backscattered  power  received  by  the  Coose  Bay  radar  is 
expected  to  be  inversely  proportional  to  the  square  of  the  distance  for  F- 
region  scatter.  For  E-region  scatter,  the  irregularities  probably  do  not  fill 
the  beam  and  the  backscattered  power  is  expected  to  be  more  closely 
proportional  to  the  cube  of  the  distance.  It  must  be  noted,  however,  that 
these  relations  are  complicated  by  a  number  of  additional  factors,  most 
notably  the  vertical  beam  pattern  of  the  antenna  array.  Figure  13  is  a 
scatter  plot  of  the  backscattered  power  (i.e.  fitted  power  from  the  ACFs)  vs 
range.  It  is  clear  that  the  density  of  points  shows  a  dependence  on  range  and 
that  the  average  power  decreases  with  increasing  range.  For  reference 
purposes  we  have  also  plotted  an  l/rJ  curve  for  the  range  300  -600  km  (where 
we  expect  the  scatter  to  come  from  the  E-region)  and  a  1/r^  curve  for  the 
range  1000  -  2100  km  (where  we  expect  the  scatter  to  come  from  the  F- 
region).  The  distribution  of  points  in  both  regions  is  reasonably  consistent 
with  the  expected  dependence  on  distance.  The  region  between  600  and  1000  km 
shows  an  unexpected  excess  of  points  with  high  power,  but  there  are  more 
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points  at  all  powers  for  this  region  and  it  contains  scatter  from  both  the  E- 
region  and  the  F-region,  which  makes  it  difficult  to  interpret  the  density  of 
points  in  this  region.  The  excess  of  points  at  high  power  for  the  closest 
ranges  (300  -  400  km)  is  probably  due  to  the  effects  of  the  antenna  beam 
pattern. 

We  now  turn  to  the  distributions  involving  the  Doppler  velocity. 
Figure  14a  is  a  scatter  plot  of  the  velocity  vs  range  for  all  the  ACFs  with 
errors  below  250  m/s  (and  power  >  0  dB).  The  most  obvious  feature  is  the 
enhanced  density  of  points  at  0  velocity.  The  most  likely  cause  for  this 
feature  is  weak  groundscatter  which,  because  the  power  was  low,  had  larger 
than  normal  spectral  widths  and  therefore  were  not  identified  as  ground¬ 
scatter.  Figure  14b  shows  the  distribution  for  only  the  higher  power  ACFs 
(.SNR  >  10  dB).  The  enhanced  density  at  0  velocity  is  no  longer  present,  but 
the  rest  of  the  distribution  is  qualitatively  the  same  as  Figure  13a,  which 
strongly  supports  the  speculation  that  the  enhancement  at  0  velocity  was  due 
to  weak  ground  scatter. 

In  Figure  15a, b  we  again  show  the  distribution  of  velocity  vs  range 
but  now  separated  in  terms  of  magnetic  activity.  The  distribution  for 
magnetically  quiet  times  (Kp  <  2-)  clearly  shows  fewer  high  velocities  than 
the  distribution  for  magnetically  active  times  (Kp  >  *  4-).  This  is 
particularly  true  at  the  closer  ranges  and  indicates  that  this  Kp  dependence 
is  significantly  more  important  in  the  E-region  than  in  the  F-region. 

Since  the  velocity  distribution  depends  on  the  high-latitude 
convection  patterns  we  expect  the  range  vs  velocity  plots  to  vary  with  local 
time.  Before  examining  this  relationship,  however,  we  must  first  consider  how 
the  geographic  orientation  of  the  Goose  Bay  radar  affects  the  observations. 
The  bore-sight  of  the  radar  points  approximately  5  degrees  east  of  geographic 
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north  rather  than  toward  the  invariant  pole  or  the  magnetic  pole.  It  is  not 
entirely  clear  which  pole  best  describes  the  plasma  convection  pattern  at  high 
latitudes.  If  we  use  the  invariant  pole  (78.5”  N,  291.0"  E)  the  direction  of 
the  pole  is  approximately  9®  W  of  the  radar  bore-sight  direction.  If  we  use 
an  eccentric  dipole  (81.7®  N,  271.8®  E)  the  direction  of  the  pole  is  about  13® 
W  and  if  we  use  the  magnetic  dipole  (76®  N,  259®  E)  the  direction  of  the  pole 
is  38.6*  W. 

Figure  16a  shows  one  model  of  the  high  latitude  convection  pattern 
appropriate  for  a  southward  pointing  interplanetary  magnetic  field  (Heelis,  et 
al.  1980).  Superimposed  on  the  convection  pattern  is  the  radar  field  of  view 
at  0,  6,  12,  and  18  UT.  In  this  view,  the  angle  at  the  radar  between  the 
geographic  and  geomagnetic  poles  is  taken  to  be  approximately  13®.  In  the 
dusk  time  sector  the  plasma  flow  is  primiarily  perpendicular  to  the  radar  line 
of  sight,  but  with  a  component  that  is  toward  the  radar  for  most  azimuths. 
Only  the  western-most  azimuth  will  see  a  negative  component.  Thus  there 
should  be  a  bias  toward  positive  velocities  for  the  18-24  UT  period.  In  the 
dawn  sector  the  situation  is  reversed.  In  the  noon  and  midnight  sectors  the 
rotation  of  the  plasma  convection  from  sunward  to  anti-sunward  (and  vice- 
versa)  should  produce  a  strong  bias  toward  negative  velocities  in  the  noon 
sector  or  a  bias  toward  positive  velocities  for  the  midnight  sector.  In  both 
cases,  the  rotation  of  the  radar  with  respect  to  the  magnetic  pole  reduces  the 
bias  from  what  would  be  expected  in  the  symmetric  case. 

Figures  17a-d  show  the  range  vs  velocity  distributions  for  the  four 
time  sectors.  The  dawn  and  dusk  distributions  (17b, d)  show  the  expected 
negative  and  positive  velocity  biases.  Figure  17a,  the  distribution  for  the 
midnight  sector,  shows  a  relatively  symmetric  pattern  with  no  evidence  of  a 
bias  toward  positive  velocities.  The  orientation  of  the  radar  would  tend  to 
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decrease  the  positive  bias  for  material  coming  out  of  the  polar  cap,  but  that 
is  insufficient  to  explain  these  observations.  One  possible  explanation  for 
this  discrepancy  is  that  the  convection  pattern  shown  in  Figure  16a  does  not 
include  the  effect  of  the  Harang  discontinuity  nor  does  it  include  the  effect 
of  the  interplanetary  magnetic  field  on  the  size  and  orientation  of  the 
convection  cells.  Figure  16b  shows  a  pattern  (Heelis,  1980)  which  may  be  more 
typical.  Some  convection  models  (e.g.  Heppner,  1984;  Friis-Christensen,  et 
al . ,  1985)  show  a  consistent  clockwise  shift  (for  all  values  of  By)  of  the 
basic  convection  pattern.  As  can  be  seen  in  Figure  16b,  this  has  little 
effect  on  what  the  radar  would  observe  in  the  dawn  and  dusk  sectors,  but  could 
strongly  affect  what  is  observed  in  the  noon  and  midnight  sectors.  The 
pattern  in  Figure  16b  would  tend  to  show  a  much  more  symmetrical  Doppler 
velocity  pattern  around  midnight.  In  addition,  By  effects  would  tend  to  shift 
the  pattern  around  the  pole.  The  combined  effects  may  be  sufficient  to  reduce 
the  expected  positive  velocity  bias  to  the  point  where  the  bias  is  no  longer 
observable.  Figure  17c  shows  the  expected  bias  toward  negative  velocities  but 
there  are  fewer  events  represented  in  this  figure  compared  to  17a.  Thus, 
although  the  velocity  pattern  is  consistent  with  either  Figures  16a  or  16b,  it 
cannot  be  taken  as  definitive  support  for  either  model. 

Finally,  we  consider  the  type  of  velocity  spectra  observed.  For  the 
purpose  of  this  study  the  parameter  which  characterizes  a  spectrum  is  the 
spectral  width.  As  was  pointed  out  in  the  discussion  of  Figure  12,  there  is 
an  inverse  relationship  between  the  backscattered  power  and  the  spectral 
width.  This  relationship,  however,  is  primarily  due  to  the  increased 
influence  of  noise  at  low  powers.  There  remains  a  question  whether  there  is  a 
real  physical  relationship  between  the  two  quantities.  Since  it  is  likely 
that  different  instability  mechanisms  operate  in  the  E-  and  F-regions,  any 
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physical  dependence  of  Che  spectral  width  on  the  backscattered  power  (and 
therefore  on  Che  fluctuation  level  of  the  ionospheric  irregularities)  might  be 
expected  to  be  different.  We  therefore  examine  the  distribution  of  spectral 
width  vs  power  for  the  E-region  and  F-region  separately  (Figures  18a, b)  and, 
as  can  be  seen,  the  two  distributions  are  virtually  identical. 

It  is  well  known  that  at  VUF  frequencies,  the  width  of  the  spectra 
tend  to  be  larger  when  the  Doppler  velocity  is  nearly  zero.  We  therefore  look 
at  the  distribution  of  spectral  width  vs  Doppler  velocity.  Figures  19a, b  show 
this  distribution  for  the  E-region  and  F-region.  The  most  obvious  feature  is 
the  two  concentrations  of  points  around  ±  450  m/ s  in  the  E-region  distribution 
(Figure  19a).  These  clusters  of  points  also  show  a  limited  range  of  spectral 
widths  from  about  30  m/s  to  125m/s.  If  these  two  E-region  clusters  are 
ignored  the  E-  and  F-region  distributions  are  quite  similar.  In  particular, 
che  E-region  distribution  does  not  show  an  enhancement  in  comparison  to  the  F- 
region  in  the  number  of  spectra  with  large  widths  at  nearly  zero  velocity. 
Both  distributions  show  a  greater  number  of  spectra  with  large  widths 
occurring  at  low  velocities,  but  this  is  probably  due  to  the  fact  that  there 
are  a  greater  number  of  points  total  at  low  velocities. 

We  next  examine  the  spectral  width  as  a  function  of  magnetic 
activity  (Kp).  Figure  20a, b  show  the  distribution  of  points  in  spectral  width 
and  range  for  quiet  conditions  (a)  and  disturbed  conditions  (b).  In  order  to 
limit  the  problem  of  noise  broadening  the  spectra,  we  have  plocted  only  those 
points  with  power  above  6  dB.  At  the  higher  ranges  (F-region)  the 
distributions  for  quiet  and  disturbed  conditions  are  very  similar,  but  at  che 
low  ranges  (E-region)  they  are  radically  different.  The  E-region  spectra 
under  quiet  conditions  show  very  narrow  Doppler  spectra,  with  most  of  the 
points  having  spectral  widths  of  less  than  60  m/s.  Conversely,  under 
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disturbed  conditions  virtually  all  the  spectra  have  widths  in  excess  of  60m/s. 

To  get  a  clearer  picture  of  the  spectral  width  distributions  we  have 
made  histograms  of  the  frequency  of  occurrence  of  spectral  width  for  the  £- 
and  F-regions  (300-600  km,  900-2100  km)  for  quiet  and  disturbed  periods. 
These  distributions  are  shown  in  Figures  21a-d.  As  can  be  clearly  seen,  the 
E-region  distributions  are  nearly  complementary.  A  comparison  of  the  E-region 
and  F-region  distributions  for  disturbed  conditions  shows  a  marked  similarity 
between  the  two,  although  on  average  the  F-region  spectra  are  somewhat  wider 
than  the  E-region  spectra.  In  comparing  the  F-region  distributions  for  quiet 
and  disturbed  conditions,  we  note  that  the  quiet  period  distribution  appears 
to  consist  of  two  components,  a  narrow  component  similar  to  the  E-region 
distribution  plus  a  broad  spectrum  component  similar  to  the  disturbed  period 
distribution.  One  might  speculate  that  in  the  F-region  there  are  two  differ¬ 
ent  irregularity  producing  processes,  one  of  which  operates  under  all  magnetic 
conditions  and  another  which  operates  only  under  quiet  conditions. 

Another  question  that  must  be  considered  is  whether  or  not  the 
spectral  widths  vary  as  a  function  of  local  time.  Since  the  results  of  the 
width  variation  with  Kp  showed  distinct  differences  between  the  E-region  and 
F-region  behaviors,  we  will  continue  to  separate  the  distributions  by 
region.  Figure  22a-d  show  the  spectral  width  distributions  for  the  4  time 
sectors:  midnight,  dawn,  noon,  and  dusk,  in  the  E-region.  There  is  no  marked 
difference  in  these  distributions.  The  12-18  UT  (noon)  sector  looks  different, 
but  there  are  very  few  points  and  they  come  primarily  from  a  single  event  so 
the  apparent  difference  is  not  statistically  meaningful.  Figure  23a-d  show 
the  spectral  width  distributions  for  the  4  time  sectors  in  the  F-region.  Here 
there  are  very  clear  differences.  In  particular,  the  narrow  spectra  clearly 
come  preferentially  from  the  dusk  and  midnight  sectors.  In  particular  the 
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width  distribution  in  the  dusk  sector  (18-24  UT)  shows  a  distinct  absence  of 
wide  spectra  compared  to  the  other  times.  Combined  with  the  results  in 
Figures  21c, d  this  suggests  that  the  narrow  spectra  are  associated  with  the 
afternoon  auroral  zone  F-region  under  quiet  conditions.  To  further 
investigate  this  supposition  we  examine  the  F-region  distributions  for  the  18- 
24  UT  period  separated  into  quiet  and  disturbed  magnetic  conditions.  (Figure 
24).  The  results  are  clear.  The  quiet  period  spectra  are  completely 
dominated  by  the  narrow  spectra  while  the  disturbed  period  distribution  shows 
very  few  narrow  spectra.  Indeed,  the  latter  distributions  are  very  similar  to 
the  results  seen  in  the  E-region  distributions  in  Figure  21a, b. 

5.0  DISCUSSION 

5.1  High  Latitude  Convection  Pattern 

As  was  pointed  out  in  section  4.3,  the  convection  pattern  in  the 
dawn  and  dusk  sectors  are  not  highly  sensitive  to  the  IMF.  However,  the 
pattern  at  noon  and  midnight  can  depend  very  strongly  on  both  the  B  and  B 
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components.  For  this  reason  a  statistical  survey,  such  as  the  one  reported 
here,  is  much  more  likely  to  give  coherent  and  consistent  results  for  the  dawn 
and  dusk  sectors  than  for  the  noon  and  midnight  sectors.  A  great  deal  of 
modeling  effort  has  been  made  to  determine  how  the  IMF  affects  the  convection 
pattern  (e.g.  Heppner,  1977;  Potemra  et  al.,  1979;  Heelis  et  al.,  1980; 
Potemra  et  al.,  1984;  Zanetti  et  al.,  1984;  Heelis,  1984;  Reiff  and  Burch, 
1985;  Friis-Christensen,  1985).  The  B^  component  has  two  primary  effects  on 
Che  convection,  a  rotation  about  the  magnetic  pole  and  a  change  in  the 
relative  sizes  of  the  two  convection  cells  (e.g.  Heelis  et  al.,  1980,  Potemra 
et  al.,  1979,  Reiff  and  Burch,  1985).  For  some  models  Che  convection  pattern 
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for  negative  By  is  the  mirror  image  of  the  pattern  for  positive  By  (Burch  et 
al.  1985$  Reiff  and  Burch,  1985;  Potemra  et  al.,  1984;  Zanetti  et  al . , 
1984).  For  others  (Heppner,  1977$  Friis-Christensen  et  al.,  1985)  this  is  not 
true.  In  addition,  some  models  include  a  dependence  on  Bx  (Heelis,  1984; 
Reiff  and  Burch,  1985).  Recently,  extensive  efforts  have  been  made  using  the 
Sondrestrom  radar  to  observe  the  effects  of  the  IMF  on  the  high  latitude 
convection  pattern  (de  la  Beaujardiere  et  al.,  1985$  de  la  Beaujardiere  et 
al.,  1986).  They  found  that  "the  IMF  By  component  strongly  influences  the 
night  side  polar  convection"  and  that  none  of  the  models  offer  a  very  good 
match  to  the  observations. 

An  additional  complication  results  from  what  are  known  as  UT 
effects.  The  precise  details  of  the  plasma  convection  depend  not  only  on  the 
magnetospheric  processes  which  drive  the  convection,  but  on  the  ionospheric 
response  as  well.  An  important  factor  in  the  ionospheric  response  is  the 
conductivity  (both  Hall  and  Pedersen),  which  is  strongly  influenced  by  the 
solar  zenith  angle  and  hence  depend  on  season  and  local  solar  time  (or  UT) 
rather  chan  the  local  magnetic  time.  Thus,  the  convection  pattern  seen  at 
different  geographic  locations  may  show  significant  differences  even  when 
comparing  identical  magnetic  local  times  and  magnetic  conditions. 

The  Goose  Bay  radar  cannot  usually  see  far  enough  into  the  polar  cap 
to  determine  details  about  the  convection  across  the  cap.  In  particular,  this 
means  it  is  unlikely  to  be  able  to  see  the  3  or  4-cell  patterns  that  have  been 
suggested  for  northward  Bz  conditions.  In  general  the  Goose  Bay  radar  results 
are  compatible  with  the  standard  two-cell  convection  patterns.  On  a 
statistical  basis  the  results  near  noon  and  midnight  show  less  of  the  anti- 
sunward  flow  than  might  be  expected,  but  when  the  complications  of  By  effects 
and  UT  effects  are  included  the  results  are  not  inconsistent  with  the  various 
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convection  models.  Also  specific  event  studies  often  yield  Doppler 
measurements  consistent  with  antisunward  flow.  A  more  detailed  study  wll 
have  to  be  made  in  the  future  in  which  the  Coose  Bay  data  can  be  separated  by 
IMF  conditions. 

5.2  Spectral  Characteristics 

The  observed  Doppler  spectra  have  a  wide  distribution  of  spectral 
widths.  These  widths  range  from  nearly  zero  (i.e.  no  decorrelation  of  the  ACF 
over  the  50  ms  multipulse  duration)  to  over  1000  m/s.  Overall  there  are  no 
dramatic  differences  between  the  E-region  and  F-region  spectra,  but  there  are 
a  few  interesting  differences.  The  E-region  spectra  appear  to  have  a  small 
additional  population  of  spectra  with  velocities  near  the  ion-acoustic 
velocity  (-  -450  m/s).  These  spectra  have  widths  from  about  30  to  100  m/s. 
We  also  note  that  beyond  Doppler  velocities  of  about  ±  500  m/s  the  number  of 
E-region  spectra  decreases  markedly,  although  spectra  are  observed  with  mean 
Doppler  velocities  as  high  as  2000  m/s.  In  contrast,  the  F-region 
distribution  of  the  mean  Doppler  velocities  shows  a  a  more  gradual  decrease  in 
occurrence  at  the  higher  velocities. 

Another  interesting  class  of  spectra  are  found  in  the  dusk  sector 
during  periods  of  low  magnetic  activity.  These  spectra  have  narrow  widths 
(-  50  m/s)  and  usually  have  Doppler  velocities  in  the  range  from  -300  to  +300 
m/s  (although  in  some  cases  the  elocities  are  over  500  m/s)  They  occur  at 
ranges  from  600  to  1800  km  and  are  presumably  from  the  F-region.  Because  of 
their  similarity  to  the  narrow  spectra  often  observed  in  the  E-region  we  must 
consider  the  possibility  that  these  echoes  come  from  a  ground  to  F  region  to 
ground  to  E-region  (and  back)  propagation  mode  and  are  in  fact  E-region 
irregularities.  There  is,  however,  no  evidence  in  the  cases  studied  for  a 
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ground  scatter  signal  equatorward  of  the  ionospheric  signal.  In  addition,  it 
seems  unlikely  that  such  a  complicated  propagation  mode  would  occur  often.  In 
this  survey  we  have  found  these  narrow  spectra  in  23  separate  events,  which 
supports  the  idea  that  these  echoes  are  indeed  from  F-region  irregularities. 
The  limited  range  of  Doppler  velocities  for  these  spectra  is  probably  not 
significant,  since  these  spectra  occur  only  during  quiet  magnetic  times  when 
the  ionospheric  electric  fields  would  tend  to  be  low. 

5.3  Implications  for  QTH  Radar  Systems 

One  thing  that  has  become  clear  in  'his  statistical  survey  is  that 
nearly  every  variety  of  spectrum  and  Doppler  velocity  is  possible.  There  is  no 
single  feature  shared  by  all  ionospheric  clutter.  Both  very  narrow  and  very 
wide  spectra  are  found.  Wide  spectra  can  occur  at  high  Doppler  velocities  as 
well  as  at  very  low  Doppler  velocities.  In  the  latter  case  a  single  volume  of 
space  can  simultaneously  exhibit  significant  backscattered  power  at  both 
positive  and  negative  Doppler  shifts.  Thus  an  entire  range  of  Doppler 
velocities  from  -500  to  +500  m/s  can  be  obscured  by  ionsopheric  clutter.  One 
cannot,  however,  avoid  this  problem  simply  by  looking  for  targets  with  very 
narrow  spectra,  since  ionospheric  clutter  with  narrow  spectra  is  also  quite 
cooinon  • 

Ionospheric  clutter  is  particularly  common  during  the  late  evening 
and  night  hours.  Nevertheless,  clutter  may  appear  at  any  time.  The  best  time 
for  observing  ground  scatter  is  the  period  around  local  noon.  However,  even 
in  the  noon  sector,  over  802  of  the  scatter  events  we  have  observed  were 
classified  as  ionospheric  clutter  or  ionospheric  clutter  plus  ground 
scatter.  The  probability  of  observing  ionsopheric  clutter  is  also  dependent 
on  magnetic  activity  and  is  much  more  likely  to  occur  during  disturbed 
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conditions.  This  also  implies  Chat  ionospheric  clutter  will  be  more  common 
during  periods  of  high  sunspot  number.  The  present  study  only  includes  data 
taken  near  the  sunspot  minimum  of  the  eleven  year  solar  cycle.  The  amount  of 
ionospheric  clutter  will  probably  be  much  greater  near  sunspot  maximum. 


37 


6.0 


REFERENCES 


de  la  Beaujardiere,  0.,  V.  B.  Wickwar  and  J.  D.  Kelly,  Effect  of  the 

interplanetary  magnetic  field  Y  component  on  the  high-latitude  nightside 
convection,  Geophys.  Res.  Lett.,  12,  461,  1985. 

de  la  Beaujardiere,  0.,  V.  B.  Wickwar  and  J.  H.  King,  Sondrestrom  Radar 
observations  of  the  effect  of  the  IMF  B^  component  on  polar  cap 

convection,  1986. 

Elkins,  T.  J.,  A  model  for  high  frequency  radar  auroral  clutter,  RADC-TR-80- 
122,  1980. 

Farley,  D.  T.,  Multiple-pulse  incoherent-scatter  correlation  function 
measurements,  Radio  Sci.,  7,  661,  1972. 

Friis-Christensen,  E.,  Y.  Kamide,  A.  D.  Richmond,  and  S.  Matsushita, 
Interplanetary  magnetic  field  control  of  high-latitude  electric  fields  and 
currents  determined  from  Greenland  magnetometer  data,  J.  Geophys.  Res., 
90,  1325,  1985. 

Greenwald,  R.  A.,  K.  B.  Baker,  R.  A.  Hutchins,  and  C.  Hanuise,  An  HF  phased- 

array  radar  for  studying  small-scale  structure  in  the  high  latitude 

ionosphere,  Radio  Sci.,  20,  63,  1985. 

Hanuise,  C.,  R.  A.  Greenwald,  and  K.  B.  Baker,  Drift  motions  of  very  high 
latitude  F-region  irregularities,  J.  Geophys.  Res.,  90,  9717,  1985. 

Heelis,  R.  A.  and  W.,  B.  Hanson,  High-latitude  ion  convection  in  the  nighttime 
F-region,  J.  Geophys.  Res.,  85,  1995,  1980. 

Heelis,  R.  A.,  J.  C.  ft>ster,  0.  de  la  Beaujardiere  and  J.  Holt,  Multistation 
measurements  of  high-latitude  ionospheric  convection,  J.  Geophys ,  Res . , 
88,  10111,  1983. 


-  38  - 


Heppner,  J.  P.,  Empirical  models  of  high-latitude  electric  fields,  J.  Geophys. 


Res.,  82,  1115,  1977. 

Heppner,  J.  P.,  Polar-cap  electric  field  distributions  related  to  Che 
interplanetary  magnetic  field  direction,  J.  Geophys.  Res.,  77,  4877,  1972. 

Nielsen,  E.,  C.  I.  Haldoupis,  B.  G.  Fejer,  and  H.  M.  Ierkic,  Dependence  of 
auroral  power  spectra  variations  upon  electron  drift,  J.  Geophys.  Res., 
89,  253,  1984. 

Potemra,  T.  A.  ,L.  J.  Zanetti,  P.  F.  Bythrow,  and  A.  T.  Y.  Lui ,  By-dependent 
convection  patterns  during  northward  interplanetary  magnetic  field,  J. 
Geophys.  Res.,  89,  9753,  1984. 

Potemra,  T.  A.,  T.  Iijima,  and  N.  A.  Saflekos,  Large-scale  characteristics  of 
Birkeland  currents,  Dynamics  of  the  Magnetosphere,  165,  1979. 

Reiff,  P.  H.  and  J.  L.  Burch,  IMF  By-dependent  plasma  flow  and  Birkeland 
currents  in  the  dayside  magnetosphere.  2.  A  global  model  for  northward  and 
southward  IMF,  J.  Geophys.  Res.,  90,  1595,  1985. 

Walker,  A.  D.  M. ,  R.  A.  Greenwald,  and  K.  B.  Baker,  Determination  of  the 
fluctuation  level  of  ionospheric  irregularities  from  radar  backscatter 
measurements,  Radio  Sci.,  1986. 


FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


A  comparison  of  the  different  pulse  techniques.  (a)  long  pulse, 
(b)  repetitive  pulse,  (c)  variable  lag  double  pulse,  (d) 
multiple  pulse. 

A  comparison  of  the  two  types  of  least-squares  fits  to  the  power 
profile  of  the  autocorrelation  function  (ACF).  (a)  An  example  in 
which  the  Gaussian  fit  is  superior,  (b)  An  example  in  which  the 
exponential  fit  is  superior.  Note  that  (b)  is  the  same  range  as 
(a)  and  only  5  seconds  later  in  time. 

A  comparison  of  the  decorrelation  parameter  \,  for  ionospheric 
backscatter  and  ground  scatter. 

The  probability  of  observing  a  backscattered  signal  3  dB  or  more 
above  the  noise  level  for  four  time  periods,  0-6  UT  (local 
midnight  sector),  6-12  UT  (local  dawn),  12-18  UT  (local  noon), 
and  18-24  (local  dusk). 

Type  of  scatter  as  a  function  of  time.  For  an  event  to  be 
"ionospheric"  at  least  25%  of  the  Doppler  spectra  (or  ACFs)  must 
be  identified  as  ionospheric.  For  an  event  to  be  ground  scatter 
at  least  25%  of  the  spectra  must  be  identified  as  ground 
scatter. 

The  type  of  scatter  as  a  function  of  magnetic  activity  (as 
defined  by  the  index).  Note  that  there  are  very  few  events 
with  high  Kp  (>  4-). 

The  probability  of  observing  ionospheric  scatter  from  the  D- 
region  and  from  the  F-region  as  a  function  of  time.  An  event  is 
E-region  scatter  if  at  least  25%  of  the  Doppler  spectra  came 
from  ranges  less  than  600  km  from  the  radar.  An  F-region  event 
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Figure  8. 


Figure  9. 


Figure  10. 


Figure  11. 


Figure  12. 


Figure  13. 


Figure  14. 


is  one  in  which  at  least  25Z  of  the  spectra  come  from  ranges 
greater  than  900  km.  Events  where  the  scatter  erne  primarily 
between  600  and  900  km  could  not  be  categorized. 

Probability  of  observing  ionospheric  scatter  from  the  E-region 
and  F-region  as  a  function  of  Kp. 

Scatter  plot  of  the  fitted  power  (C)  from  the  least-squares  fit 
to  the  ACF  vs.  the  lag-o  power  for  the  ACF. 

Scatter  plot  of  the  estimated  error  in  the  mean  Doppler  velocity 
vs.  (a)  the  lag-o  power  of  the  ACF  and  (b)  the  fitted  power. 
Scatter  plot  of  the  estimated  error  in  the  mean  Doppler  velocity 
vs.  the  spectral  width  (in  m/s)  as  determined  from  the 

decorrelation  parameter,  X. 

Scatter  plot  of  the  spectral  width  vs.  the  fitted  power  in  the 
ACF.  The  sudden  decrease  in  the  density  of  points  at  low  power 
and  low  width  is  due  to  the  removal  of  spectra  which  were 
identified  as  ground  scatter  signals. 

Fitted  power  vs.  range  (in  km)  from  the  radar.  For  echoes  from 
the  thin  E-region,  the  power  should  be  inversely  proportional  to 
the  cube  of  the  distance.  For  the  F-region,  the  power  should  be 
inversely  proportional  to  the  square  of  the  distance. 

Scatter  plot  of  the  range  (in  km)  vs.  the  mean  Doppler 
velocity.  (a)  Scatter  plot  for  all  points  with  velocity  error 
less  than  250  m/s  and  signql-to-noise  ration  >  0  dB.  (b) 
Scatter  plot  for  only  the  points  with  velocity  error  <  250  m/s 


and  SNR  >  10  dB. 
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Figure  IS. 


Figure  16. 


Figure  17. 


Figure  18. 


Figure  19. 


Figure  20. 


Figure  21. 


Scatter  plot  of  the  Range  vs.  the  mean  Doppler  velocity  for  (a) 
All  points  with  velocity  error  <  250  m/s  and  Kp  <  2-  (magnetic 
quiet  times)  and  (b)  All  points  with  veLocity  error  <  250  m/s 
and  Kp  >  4-  (magnetically  disturbed  times). 

Examples  of  proposed  convection  patterns  (Heelis  and  Hanson, 

1980)  in  relation  to  the  f ield-of-view  of  the  Goose  Bay  radar. 

(a)  Symmetric  pattern  for  B„  <  0  and  B„  =  0.  (b)  Asymmetric 
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pattern  for  Bz  <  0  and  B^  >  0.  Mote  that  some  convection  models 
show  a  similar  pattern  independent  of  the  sign  of  B^. 

Range  vs.  velocity  for  all  points  with  velocity  error  <  250  m/ s 
and  time  in  the  (a)  0-6  UT  range,  (b)  6-12  UT  range,  (c)  12-18 
UT  range,  and  (d)  18-24  UT  range. 

Spectral  width  vs.  the  fitted  power  for  (a)  The  E-region  (range 
300-600  km)  and  (b)  The  F-region  (range  900-2100  km). 

Spectral  width  vs.  the  mean  Doppler  velocity  for  (a)  The  E- 
region  (300-600  km)  and  9b)  the  F-region  (900-2100  km). 

Spectral  width  vs.  the  Range  for  all  spectra  with  signal-to- 
noise  ration  >  6  dB  (to  eliminate  spectra  which  are  broader  than 
normal  because  the  signal  is  weak)  and  velocity  error  <  250  m/s 
and  (a)  Kp  <  2  (quiet  periods)  and  (b)  Kp  >  4-  (disturbed 
periods). 

-^Spectral  width  distribution  for  all  points  with  signal-to-noise 
ration  >  6  dB  and  velocity  error  <  250  ml s  and  (a)  E-region 
echoes  (300-600  km)  under  quiet  (Kp  <  2-)  conditions,  (b)  E- 
region  echoes  under  disturbed  (Kp  >  4-)  conditions,  (c)  F-region 
echoes  (900-2100  km)  under  quiet  conditions,  (d)  F-region  echoes 
under  disturbed  conditions. 
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Figure  22. 


E-region  spectral  width  distribution  as  a  function  of  time,  (a) 
0-6  UT  (local  midnight),  (b)  6-12  (local  dawn),  (c)  12-18  UT 
(local  noon),  and  (d)  18-24  UT  (local  dusk.). 

Figure  23.  F-region  spectral  width  distribution  as  a  function  of  time,  (a) 
0-6  UT,  (b)  6-12  UT,  (c)  12-18  UT,  and  (d)  18-24  UT. 
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APPENDIX 


The  following  is  a  brief  description  of  Che  mechod  Co  generate  Che  4 
and  6  dimensional  contingency  Cables  described  in  sections  4.2  and  4.3. 

In  doing  Che  staciscicaL  analysis  of  Che  radar  data,  it  is  necessary 
to  write  a  computer  program  Chat  will  put  n  variables  of  Che  radar 
observations  such  as  time,  range,  power,  etc.,  into  an  n-dimensional  histogram 
or  array  and  then  project  this  array  to  (n-l)-dimensional  arrays  and  project 
the  resulting  arrays  to  (n-2)-dimensional  arrays.  This  process  will  continue 
until  everything  is  reduced  to  1-dimensional  array. 

The  program  would  be  very  simple  and  straight  forward  to  write,  if 
there  were  only  3  or  fewer  variables,  thus  using  only  arrays  of  3-dimension  or 
smaller.  However,  the  program,  if  not  carefully  planned,  could  explode  into 
pages  and  pages  of  do  loops  and  write  statements  as  more  variables  were 
added.  This  is  because  the  number  of  arrays  increases  according  to  the 
recursive  function  numarray(n)  *  numarray(n-l)  x  2  ♦  1,  and  numarray(l)  =  1, 
where  numarray  is  the  number  of  arrays  produced  for  an  n-dimensional  array. 
For  example,  in  doing  4-dimensional  array  of  time,  region  type,  scatter  type, 
and  Kp  (a  magnetic  disturbance  index),  we  need  the  total  of  the  following  15 
arrays : 

1  4-dimensional  array  of  the  above  4  variables 

4  3-dimensional  arrays  (permutations  of  3  out  of  the  4  variables) 

6  2-dimensional  arrays(permutations  of  2  out  of  the  4  variables) 

4  1-dimensional  arrays 

So,  for  5-dimensional  array,  according  to  the  formula  there  are  2(15)+1  =  31 
arrays  and  for  6-dimensional  array  there  are  2(31 )+l  =  63  arrays. 

Obviously  it  would  be  impractical  to  create  63  separate  arrays  for 
each  of  the  permutations,  sum  them  up  and  print  them  out.  A  clever  algorithm 
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needs  to  be  found  for  the  program  so  that  it  is  efficient,  robust,  concise  and 
yet  easy  to  change  should  there  be  a  need  to  add  or  subtract  variables/ 
dimension.  The  following  describes  such  algorithm. 

The  above  4-dimensional  array  of  time,  region  type,  scatter  type, 
and  Kp,  will  serve  as  an  example.  These  variables  are  called  t,  rt,  st,  and 
Kp  respectively.  Furthermore,  the  analysis  calls  for  the  time  to  be  divided 
into  4  divisions,  region  type  into  4,  scatter  type  into  3,  and  Kp  into  10. 
Thus,  the  histogram  array  could  be  declared  in  fortran  as  histo(4 ,4 ,3 , 10) . 
The  simple  method  in  getting  all  the  15  arrays  would  generate  a  fortran 
program  that  looks  like  the  following: 
do  t  *  1,4 

do  rt  ■  1,4 

do  st  ■  1,3 
do  kp  -  1,10 

c  begins  listing  the  15  histograms/arrays  from 

c  the  fewest  dimensions  to  the  most  dimension 

sum_t(t)  *  sum__t(t)  +  histo(t,rt,st,kp) 
sum_rt(t)  =  sum_rt(rt)  +  histo(t,rt,st,kp) 
sum_t_st(t,st)  *  sum_t__st(t,st)  +  histo(t,rt,st,kp) 
sum_t__rt_kp(t,rt,kp)  *  sum_t__rt__kp(t ,rt,kp)  +  histo(t ,rt,st,kp) 
<etc> 
end  do 
end  do 

To  print  all  the  15  histograms,  the  above  codes  will  be  repeated 
roughly  15  times  and  therefore  will  not  be  listed  here. 


The  better  alternative  is  to  have  one  huge  array  which  contains  all 
the  15  arrays  instead  of  15  smaller  ones.  The  idea  is  to  have  just  one  array 
that  can  be  indexed  properly  such  that  it  can  be  inserted  in  a  few  do-loops  to 
give  all  the  15  arrays  and  therefore  eliminates  the  need  to  enumerate 
explicitly  each  of  the  15  individual  arrays.  (15  in  this  case,  but  it  will 
grow  to  63  for  a  6  dimensional  array  and  127  for  7  dimensional  array).  This 
big  array  is  declared  as  sum(0:4,0:4,0:10,15) .  Then, 
sumjt(t)  is  replaced  with  sum(0,0,t, 1). 
sum_rt(rt)  with  sum(0,0,rt,2). 
sumjt_st(t,st)  with  sum(0,  t , st ,  7)  etc. 

Also,  for  reasons  that  will  become  clear  later,  t  is  equivalent  to  el(l),  rt 
to  el(2),  st  to  el(3),  kp  to  el(4),  since  t  is  the  1st  element  in  the  array 
histo,  rt  is  2nd  etc.  After  doing  all  of  the  substitutions  the  listing  of  the 
above  codes  now  looks  like  the  following: 


the  first  four  do-loops  of  t,rt,st,and  kp  stay  the  same 


sum(0 ,0,el(  1 )  ,  1 )  *  sum(0 ,0 ,el(  1 ) ,  1 )  +  histo(t,rt,st,kp)  . !sura_t 

sum(0,0,el(2)  ,2)  =  sum(0,0,el(2),2)  +  histo(t,rt,st,kp)  . . .'sum  rt 


sum(0,el(l),el(3),7)  -  sum(0 ,el( 1 ) ,el*3) ,7)  +  .histo(t,rt,st,kp)  !sum_t_st 


sum(el( 1 ) ,el(2) ,el(4) , 15)  -  sum(el( 1 ) ,el( 2) ,el(4 ) , 15)  +  histo(t,rt,st,kp) 
.'sum  t  rt  kp 
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The  fourth  index  of  the  array  sum  is  numbered  from  1  to  15  for 
the  15  arrays  so  that  each  of  the  15  has  a  unique  number. 

The  next  immediate  step,  which  is  the  final  step,  is  to  find  a 
way  to  generate  the  above  list  with  do-loops.  The  key  to  doing  that  is  to 
notice  the  first  3  indices  of  the  array  el  go  from  001,  002,  003,  004, 
012,  013,  all  the  way  to  123,  124  and  234.  The  first  index  is  always  less 
than  the  second,  and  the  second  from  the  third,  except  when  the  index  is 
zero.  This  pattern  can  be  easily  generated  by  3  do  loops  as  shown  in  the 
next  few  lines. 
el(0)  -  0 
compindx  *  0 
do  indxl  ■  0,4 

do  idx2  *  indxl ,4 

if (indx2.eq.0) .or. (indx2.ne. indxl))  then 
do  lndx3  *  indx2+l,4 
compindx  *  compindx  +  1 

sum(el(indxl ) ,el(indx2) ,el(indx3) .compindx)  = 
sum(el( indxl ) ,el( indx2) ,el( indx3) .compindx)  + 
his  to2 ( t , rt , s  t , kp ) 
end  do 
end  if 
end  do 
end  do 

The  algorithm  to  print  the  arrays  looks  similar  to  the  above 
algorithm.  It  has  the  same  number  of  do-loops  and  complexity.  This  algorithm 
will  execute  almost  as  fast  as  the  simple  algorithm.  The  only  factor  that  can 
probably  slow  it  down  is  the  memory  paging  since  the  array  sum  can  get  big  as 
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new  variables  are  added.  However,  once  Che  algorithm  is  understood,  little 
modification  to  the  algorithm  will  eliminate  or  reduce  the  problem,  at  the 
expense  of  the  program  being  more  difficult  to  follow.  The  complete  listing 
of  the  subroutine  output,  which  sums  all  the  histograms  and  outputs  them,  is 
given  below. 
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1 3-Mo r-1 986  15:16:48  VAX-11  FORTRAN  V3.2-37 
1 3-Mo r-1 986  15:15:26  SYSSSYSUSER: [WING.HFR.ST/ 


0081 

subroutine  output2 

0002 

c 

0003 

c 

sun  and  output  the  4-d  histogram 

0004 

c 

the  histogrom  is  located  in  the  common  block  outdot 

0005 

c 

0006 

imp) ici t  none 

0007 

include  ‘outdot . par/I i st ' 

0008 

1 

c 

0009 

1 

c 

common  for  the  output  buffer  and  histogrom 

0010 

1 

c 

0011 

1 

integer*2  hi stol (10,4.3.8,3,4)  I6d  histogram  containing 

0012 

1 

c 

kp.w.p.v, r.t 

0013 

1 

integer*2  hi sto2(4.4,3, 10)  1  4d  histogram 

0014 

1 

c 

t , rt ,st , kp 

0015 

1 

integer*2  h i stogrnd(4,3, 10)  i  ground  scatter  histogram 

0016 

1 

c 

t.rt.kp 

0017 

1 

equivalence  (h i stol . h i sto2; 

0018 

1 

common/outdot/hi stol . hi stogrnd 

0019 

lnteger«2  t.rt.st.kp  1  indices  of  histo 

0020 

integer«2  ( ndxl , i ndx2 . i ndx3  1  indices  to  array  el 

0021 

integer*2  el (0:4)/0,0.0,0,0/  I  indices  to  sum 

0022 

equivalence  (el(1).t),(el(2),rt1,(el(3),st),(el(4).kp) 

0023 

integer*2  sum(0:4, 0:4,0: 10,0: 15/  lthe  sun  of  oil  15  histograms 

0024 

integer«2  compindx  1  the  4th  Index  of  sum 

0025 

integer*2  nel(0:4)/1 ,4.4.3,10/  i  dimensions  of  hi sto(4,4.3. 10) 

0026 

i nteger»2  i . j . k. 1 ,m  1 

0027 

character*5  coll(0:4)/‘  ‘ . ‘TIME’ . ‘REGIO’ . ‘SCATY’ . ‘KP  '/  1  label 

0028 

character*3  tob 1 e(0: 4)/ ‘  ‘ ‘ t— ' . * rt-' , ’ st-’ , ’ kp-'/  1  label  for 

0029 

characters  col(0:10,4)  1  more  label 

0030 

0031 

c 

. 

0032 

c 

initialize  the  1 abe 1 s 

0033 

c 

0034 

c 

0035 

c 

used  when  index  -  0  in  . array  sum 

0036 

c 

0037 

co 1 (0 . 1 )  <■  ‘  ‘ 

0038 

col (0.2)  -  ‘  ‘ 

0039 

col  (0.3)  >■-  ‘  ‘ 

0040 

col (0.4)  -  ’  ‘ 

0041 

c 

0042 

c 

label  for  time 

0043 

c 

0044 

col(1 .11  -  ‘0—6’ 

0045 

col (2.1 )  -  ‘ 6-1 2 ‘ 

0046 

co 1 (3 , 1 )  -  ' 12-18' 

0047 

col(4,1)  -  '18-24' 

0048 

c 

0049 

c 

label  for  region  type 

0050 

c 

0051 

col(1. 2)  -  'E' 

0052 

co 1 (2 ,2 )  -  'F' 

0053 

col (3.2)  -  'E  k  F' 

0054 

col (4,2)  -  'U' 

0055 

c 

0056 

c 

label  for  scatter  type 

0057 

c 
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output 


OUTPUT 2 


1 3-Mo r-1 986  15:16:48 
1 J-Mor-1986  15:15:26 


VAX-11  FORTRAN  V3.2-37 
SYSJSYSUSER : [ W I NG . HFR . ST 


8858 

col 

(1.3)  - 

'  I  OIOS  ■ 

8859 

col 

2.3  - 

’ GRND ‘ 

8868 

col 

(3.3)  - 

‘ION  Jc  GRND 

8061 

8862 

8863 


8874 

0875 

8876 

8877 
8078 
8879 
8888 
8081 
8082 

8883 

8884 

8885 
0886 
0087 
0088- 
8089 
8898 

8891 

8892 

8893 

8894 
0895 
0096 
0097 

8098 

8099 

8100 
0101 
0102 
0103 
0104 
0105 
8106 
0107 
0108 
0109 
0110 
0111 
0112 
0113 
8114 


label  for  kp 


8864 

col 

(1.4)  - 

•0f 

8865 

col 

(2.4)  - 

•1* 

0866 

col 

(3.4)  - 

*2' 

0867 

col 

4,4)  - 

’3' 

8868 

col 

(5.4)  - 

•41 

8869 

col 

6,4)  — 

•5’ 

8878 

col 

7,4)  - 

•6' 

0871 

col 

8.4)  - 

'7' 

0072 

col 

9.4)  - 

•8* 

8073 

col 

[10. 4)  - 

•u 

•um  I s  initialized  to  0  by  default  to  there  It  no  need  to 
initialize  it  aga i n 


sum  the  histogram 

do  t  -  1,4 

do  rt  ■  1,4 
do  et  -  1 .3 
do  kp  -  1,10 
compindx  -  0 
do  indxl  -  0,4 

do  i ndx2  -  i ndxl , 4 

i f ( ( i ndx2 . eg. 0) . or . ( i ndx2. ne . i ndxl ) ) then 
do  indx3  -  indx2+1,4 

compindx-  compindx+1  I  increment  index 

now  do  the  summation 

sum(el ( indxl) ,el ( indx2) ,el ( indx3) . 
compindx)  -  sum(e I ( i ndxl ) . e I ( i ndx2) , 
e  I  (  i  ndx3)  , comp i ndx)  4-  h(eto2( 
t.rt.et.kp) 
end  do 
end  i  f 
end  do 
end  do 
end  do 
end  do 
end  do 
end  do 

print  the  sum  histogram  out 


compindx  -  0  I  initialize  compindx  to  0 
do  indxl  -  8,4 

do  i ndx2  -  I ndxl . 4 

I f (( indx2.eq.0) . or . ( i ndx2 . ne . i ndxl ) )  then 
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0UTPUT2 


13-Mor-1986  15:16:48  VAX-1 1  FORTRAN  V3 . 2-37 
13-Mor-1986  15:15:26  SYSSSYSUSER : [KING . HFR . ST/ 


8115 

8116 

8117 

8118 
8119 
8128 
8121 
8122 

8123 

8124 
ei25 
8126 

8127 

8128 
8129 

8138 

8131 

8132 

8133 
0134 

8135 

8136 
0137 
0138 

8139 

8140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 


do  Indx3  «  lndx2+1,4 

compindx  -  compindx  4  1  I  Increment  the  4th  index 
c 

c  output  what  kind  of  hlstogrom  i«  being  outputted 

c 

wr i t e (2 , 90)  table(indx1),table(indx2),table(indx3) 

90  f ormot ( ' 1 1 . '  histoarom  of’,3o3) 

do  I  —  1  ,  ne  I  (  i  ndxl  ) 
i  f ( i ndxl . ne .8) then 
e I ( i ndxl )  —  i 
end  i  f 
c 

c  output  the  label  of  the  outer  index 

c 

wr!te(2.30)  coll(i ndx1),col(e!(i ndxl), I ndxl) 

30  f ormot ( "0 ‘ , a5 , al 0) 

c 

c  output  the  colum  labels 

c 

write(2,60)  co  1 1  ( i  ndx3)  ,  (col  (  k.  i  ndx3)  , 

#  k-1 , ne I ( I ndx3) ) 

60  f ormot ( t40, a5 ,/, t20 ,<ne I ( i ndx3)>o10) 

do  J  -  1,nel(lndx2) 
i f ( i ndx2 . ne . 0)  then 
el(indx2)  -  j 
end  i  f 
c 

c  output  the  row  labels  and  the  content  of  histogram 

c 

write(2.10)coll(indx2),eol(el(indx2), 

#  indx2),(sum(el(indx1),el(indx2),k.compindx), 

#  k  -  1 , ne I ( i ndx3) ) 

10  f ormot ( 1x,a5,o10,t18,<nel(indx3)>i10) 

end  do 
end  do 
end  do 
end  i  f 
end  do 
end  do 
return 
end 
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FIGURE  17b 
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FIGURE  17d 
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FIGURE  19a 
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FIGURE  20b 
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FIGURE  21c 
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FIGURE  22c 
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FIGURE  23a 
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FIGURE  24b 


95 


i 

|  MISSION  I 

\  .  °f  h 

§  Rome  Air  Development  Center  2 

I  l 

%  RAJDC  plans  and  executes  research,  development,  test  and  selected  3 

¥  acquisition  programs  in  support  of  Command,  Control,  Communications  % 
&  and  Intelligence  (C3!)  activities.  Technical  and  engineering  support  within  <3 

•  areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 

n  ESD  elements  to  perform  effective  acquisition  of  C3 I  systems.  The  areas  ^ 
v  of  technical  competence  include  communications,  command  and  control,  Q 
igf  battle  management,  information  procasing,  surveillance  sensors,  v 

&  intelligence  data  collection  and  handling  solid  state  sciences,  <jj 

*  electromagnetics,  and  propagation,  and  electronic,  maintainability,  and  ^ 

k  compatibility.  ^ 


